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ABSTRACT 

We use HST/ACS to study the resolved stellar populations of the nearby, nearly edge-on galaxy NGC 4244 
across its outer disk surface density break. The stellar photometry allows us to study the distribution of different 
stellar populations and reach very low equivalent surface brightnesses. We find that the break occurs at the same 
radius for young, intermediate age, and old stars. The stellar density beyond the break drops sharply by a factor 
of at least 600 in 5 kpc. The break occurs at the same radius independent of height above the disk, but is 
sharpest in the midplane and nearly disappears at large heights. These results make it unlikely that truncations 
are caused by a star formation threshold alone: the threshold would have to keep the same radial position 
from less than 100 Myr to 10 Gyr ago, in spite of potential disturbances such as infall and redistribution of 
gas by internal processes. A dynamical interpretation of truncation formation is more likely such as due to 
angular momentum redistribution by bars or density waves, or heating and stripping of stars caused by the 
bombardment of dark matter sub-halos. The latter explanation is also in quantitative agreement with the small 
diffuse component we see around the galaxy. 

Subject headings: galaxies: evolution — galaxies: halos — galaxies: individual (NGC 4244) — galaxies: 
spkal — galaxies: stellar content — galaxies: structure 



1. EVITRODUCTION 



Using photographic plates of edge-on galaxies, |van der| 
|Kruit | (1979) was the first to note that some spiral galax- 
ies have truncated disks. More recent studies with CCDs of 
galaxies at a range of inclinations have found that these trun- 
cations are not the end of the stellar disk, but are better de- 
scribed as a change in the exponential length scale of the ra- 
dial light distribution ( [Pohlen et al.|2002| . These downward 
breaks in the light profile occur in the majority of galaxies 
( [Pohlen & Trujillo,,2006) . Many different models have been 
proposed to explain the truncations or breaks but they can all 
be placed into two broad categories: models that form stars 
in the observed configuration and models that transport stars 
into the observed configuration by dynamical effects. Here 
we use observations of the resolved stellar populations of the 
edge-on galaxy NGC 4244 to constrain these types of models. 

The observed projection of a galaxy can affect the de- 
tectability of breaks in the light profiles. Breaks were orig- 
inally discovered using edge-on galaxies (van der Kruit 1979; 
|van der Kruit & Searle 1981 ) and seem to be more readily 
observed in edge-on galaxies than face-on galaxies. This can 
be partly explained by the line-of-sight projection that allows 
one to probe fainter (de-projected) surface brightnesses for 
edge-on galaxies. Furthermore, to reach faint isophotes in 
face-on galaxies, azimuthal smoothing is used to extract ra- 
dial profiles. This technique, however, smooths out the non- 
axisymmetric features (bars, spiral arms) that could result in 
sharp features in edge-on projections. 

A number of recent studies with larger samples of galax- 
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ies have tried to bridge the projection gap and have identi- 
fied a number of characteristics of breaks. |Erwin et aLlpOOS) ) 
and Pohlen & Trujillo| (2006 1 identify three types of luminos- 
ity profiles: those with a decreasing scale length beyond the 
break, those with an increasing scale length, and those with 
no break at all. Downward breaks generally occur at 2-A disk 
scale lengths, with galaxies at the short end of that range typ- 
ically exhibiting lower central surface brightnesses and larger 
scale lengths. These correllations cause the breaks to occur 
at a similar surface brightness in different galaxies ( [Kregel & 
Ivan der Kruit 2004). The breaks occur at about the same ra- 



dius independent of the height above the midplane (Pohlen 
et al .12007 ), and appear to be present already at redshifts z 



( |Perez|20(j4t|Trujillo & Pohlen|2005| . There are some indi- 
cations that disk breaks and (gaseous) warps are related, but 
there are exceptions to the rule ( van der KruitpOO? I 



As mentioned earlier, models to explain the downward 
breaks in galaxy disks come in two broad categories. In the 
first group the breaks are created by stars forming intrinsically 
in the observed distribution, either as a result of a limit in the 
gas distribution (van der Kruit|1987l l , or s ome kind of thresh- 
old suppressing star formation (e.g.,|Fal l & Efstathiou 198()| 
Kennicutt| ri989; Schaye 2004; Elmegre eir& Hunter||2006f ! 



In the second category are models that rearrange the stars 
after they form. This may occur by redistribution of angu- 
lar momentum due to bars and spiral arms (e.g., DebattistaJ 
et al.|2006 Foyle et al., in prep.), or by stripping of stars due 



to rep eated tidal interaction with other galaxies (e.g., |Gnedin| 
2003| l or dark matter sub-halos (Kazantzidis et al., in prep.). 



To constrain the large number of models in the literature 
we investigate the resolved stellar population along the ma- 
jor axis of NGC 4244 to study the age distribution of stars 
across its break. NGC 4244 is a late-type S Acd spiral galaxy 
with a maximum rotation speed of about 95 km s"' . We use a 
distance modulus of 28.20 for NGC 4244 as determined from 



the tip of the Red Giant Branch (RGB) method by |Seth et aL 
( |2005a| . The observations and data reduction are described 
in ^ the stellar distributions are investigated in ^ and we 
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Fig. 1 . — Color-magnitude diagram of tlie ACS field straddling the North- 
East break of NGC 4244. The colored polygons identify the marked popula- 
tions. Colored points show the selected stars, limited to mpgMW < 25.9 mag. 



finish with discussion and conclusions in ^ 

2. OBSERVATIONS AND DATA REDUCTION 

Our stellar photometry was derived from HST observations 
in the F606W and the F814W bands with the Advanced Cam- 
era for Surveys (ACS) Wide Field Camera. We obtained 
SNAP observations with typical total exposure times of about 
700 s in each filter split across two exposures. This data set is 
part of the GHOSTS survey of 14 galaxies with similar data 
as presented here (jde Jong et al.|2007)l. 

We used DOLPHOT (Dolphin 2000) to obtain PSF-fitting 
stellar photometry of our fields. The output of DOLPHOT 
was filtered to select only real stars using a combination of 
criteria on signal-to-noise (S/N), magnitude, and DOLPHOT 
flag and sharpness parameters. These criteria remove most 
non-stellar objects (galaxies and image defects) while leav- 
ing a large fraction (but certainly not all) stellar objects in 
the sample. Contamination from unresolved galaxies occurs 
only at low S/N, where a sharpness criterion is too poorly de- 
fined to separate all galaxies from stars. At the bright end 
there may be some Galactic foreground contamination, but 
this is expected to be less than 10 stars per 3'x3' ACS field 
for otf8i4w=23-26 mag ( Girardi et al.|2005 1. 



We ran DOLPHOT artificial star tests with at least 1 mil- 
lion stars per ACS pointing to calculate incompleteness cor- 
rections as a function of stellar magnitude, color, and crowd- 
ing. Each observed star was given an observing probability 
based on these artificial star tests, which was used to correct 
the star count radial profiles described in the next section. The 
average correction is about 2.5 in the most crowded regions, 
reducing to less than 1.1 in the outer regions. Full details of 
the GHOSTS data pipeline will be described in a forthcoming 
paper (Radburn-Smith et al., in prep.). 

3. STELLAR DENSITY PROFILES 

In Fig.[T]we show the color-magnitude diagram for the ACS 
field at the North-East end of the disk that covers the disk 
break. Several different stellar evolutionary stages can be 
identified in this diagram. We have used polygon shapes to 
identify our separation in broad evolutionary stages, but one 
should realize that some overlap in age will exist in the differ- 
ent polygons for certain ages and metallicities (see also Seth 
let al.|2005b| l. 



The magenta polygon indicates the region dominated by 
stars less than 100 Myr old. The main sequence stars are on 
the left. Barely separated from the main sequence at about 
'Mf606w-'Mf8I4w~0.1 are the blue He burning stars. The spur 
up at about mF606w-'WF8i4w^l-0 are red He burnings stars. 
These stars have been partly excluded from this bin as they 
overlap with the other populations. We refer to this bin as the 
YMS (young main sequence) population. The blue polygon 
is dominated by older main sequence and He burning stars 
of about 100-300 Myr old (labeled the MS population). The 
green polygon indicates the region dominated by Asymptotic 
Giant Branch (AGB) stars of at least 0.3 Gyr, but mostly of in- 
termediate age (1-3 Gyr), with a tail to older ages of up to 10 
Gyr The red polygon is mostly populated by metal-poor RGB 
stars. Almost all these stars will be older than 5 Gyr How- 
ever, even for a constant star formation rate population this 
polygon is dominated by ^8-12 Gyr old stars and younger 
populations can only dominate if star formation was skewed 
toward younger ages. We will use these four population boxes 
to investigate the spatial distribution of stars of different age. 

In the bottom part of Fig. |2]we show an SDSS image of 
NGC 4244 with the ACS fields overlaid. Shown in color are 
five strips used to extract star count profiles. In the top part 
of Fig. |2] we present the radial star count in these five strips 
parallel to the major axis for our four population boxes. The 
contamination subtracted star counts are presented as sym- 
bols with 2-sigma Poisson noise error bars. Downward point- 
ing arrows indicate that the Poisson noise would move these 
points below the contamination level as determined from the 
average contaminant (Galactic stars, unresolved galaxies, im- 
age defects) density in the outermost major and minor axis 
fields. 

Most conspicuous in the comparison of the different mid- 
plane profiles is that they all show a break in their exponential 
slope at about -420" (-9.0 kpc), independent of stellar popu- 
lation. Whatever causes this truncation, it affects stars of very 
different ages in the same way. The projected scale length 
in the midplane changes by at least a factor 4 for RGB and 
AGB stars, and even more so for the young stars as they have 
a much larger inner scale length. This flatter distribution for 
young populations could be caused by star formation mainly 
occurring in a ring or in wide spiral arms. The radial distri- 
bution is not symmetric in both directions — especially for 
young stars, supporting the aforementioned spiral arm inter- 
pretation — and our South- West field is not at a large enough 
radius to detect the break on that side. 

The MS and RGB populations have steep midplane surface 
density profiles beyond the break, dropping by at least a factor 
of 600 in 5 kpc. For the YMS and AGB profiles we cannot ac- 
cess the exact behavior across the break because of the much 
lower counts so that the profiles get lost in the gap between 
two ACS fields at -500" to -600". However, the upper limits 
in the outer fields show that a significant break must occur for 
these populations as well. 

The RGB star count profiles agree well with luminosity pro- 
files extracted from the SDSS image, and by lining them up, 
we find that we reach about 31 iAB-mag arcsec"^ equivalent 
surface brightness. Our measurements also agree well with 
previous inte grat ed light measure ments by van der Kruit & 
Searle| ( [T98T] l and |Fry et aLldl999| l. These authors find a trun 

^K^ii*- znr\" <-Un« /ior\" k^^^i^ 



cation radius of about 570", larger than our 420 break ra- 
dius. However, their truncation measurements are the average 
of the shorter North-East break radius observed here and the 
larger South- West break radius, and the edge-on projection of 
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Fig. 2. — Top) Incompleteness corrected star count profiles of NGC4244. 
Profiles are shown for a 70" wide strip along the midplane, and for 35"-wide 
snips offset 50" and 80" above and below the plane. The different symbols 
and colors used for the 5 strips are identified in the legend. Contamination- 
subtracted star count surface density (S) profiles are presented for the four 
CMD regions identified in Fig.^as marked, where for clarity the AGB, MS, 
and YMS populations have been offset by —2, —6, and —10 respectively in 
log(S). The error bars are 2-sigma Poisson noise uncertainties, with arrows 
indicating that Poisson uncertainty will move the counts below the back- 
ground level. Two-sigma upper limits are plotted for bins in which we detect 
fewer than 3 stars. For reference we show the equivalent SDSS ('-band surface 
brightnesses on the right axis, which was obtained by lining up the SDSS i- 
band integrated light profiles with the RGB stellar density profiles. Bottom) 
SDSS image of NGC4244 rotated by 42.7°, such that the North-East is to 
the left. Overlaid are the positions of our ACS fields in gray and the position 
of the strips used to extract the profiles above, in the same color coding. 



the truncation causes deviations from an exponential profile 
at radii smaller than the truncation radius. 

At about -640" the MS and RGB star count profiles seem 
to change scale length again to a shallower slope. This occurs 
close to our detection limit and the shape and slope of this 
outer profile is poorly constrained due to the uncertainty in the 
contamination level. This might be a diffuse halo, akin to our 
detection on the minor axis ( |Seth et al.| 2007 ), as this compo- 
nent shows only marginal concentration toward the midplane. 

Figure |2]also highlights the behavior of the break as a func- 
tion of height above the midplane. Inside the break, the pro- 



files at each of the vertical offset positions have very similar 
radial distributions for a given population. The density dif- 
ferences between the different vertical positions are larger for 
the younger populations due to their shorter scale height ( |Seth| 
|et al.|20 d5b ). The break occurs at roughly the same radius in 
the 50" offset profiles as in the midplane profiles. Beyond 
-420" the slope of the profile gets shallower because all pro- 
files beyond the break, independent of the height above the 
disk, converge to the same star count level at about 640" from 
the center At 80" above and below the disk there is no longer 
a break visible as the profile slope beyond the break has be- 
come equal to the inner slope. 

4. DISCUSSION AND CONCLUSIONS 

Our HST/ACS star count data for NGC4244 allow us to 
derive the following observations: 

i) The North-East major axis radial star count profiles show a 
sharp break with two exponential sections. The stellar surface 
density drops by at least a factor of 600 before reaching a dif- 
fuse halo or the background contamination level. 

ii) The surface density profiles on the major axis for the 
young, intermediate age, and old populations all show a break 
at the same radius. 

iii) The break occurs at the same radius independent of height 
above the midplane. The inner scale length is relatively con- 
stant, but the outer slope becomes gradually shallower with 
increasing disk height. This means that in face-on projections 
the break will be less pronounced. 

The observation that all stellar populations, regardless of 
age, have similar breaks in the midplane profiles puts the 
strongest constraint on models of break formation, assum- 
ing that the NGC 4244 break is representative of other breaks. 
Most notably, any model where breaks are caused by a sud- 
den change in star formation rate will now have to explain 
why the radius at which this occurs is stable over a very long 
time. The model will also have to explain why the breaks 
in the different populations have similar shapes, whereas one 
would expect that a very old break would be smoothed by the 
dynamical heating that we know occurs in older populations. 

For instance, a model like that of [van der Kruit| ( |1987| ), 
in which the break results from the specific angular momen- 
tum distribution of a collapsing gaseous sphere that limits the 
maximum angular momentum available to redistribute gas in 
the outer disk, has several problems. Over time, more and 
more gas will decouple from the Hubble flow with different 
distributions of specific angular momentum, and it is hard to 
conceive that the cutoff radius would never change. Further- 
more, the formation of bars and spiral arms redistributes an- 
gular momentum from the central regions to the outer disks 
so that the current specific angular momentum distribution no 
longer nece ssarily reflects the initial distribution (e.g., Debat^ 
tista et aL|20 06 ). Finally, as the Hi gas distribution of many 
galaxies is significantly more extended than their stellar dis- 
tribution, it would seem that higher angular momentum gas is 
often akeady present (although this is currently not true for 
NGC 4244; Oiling 1996). 

Models that rely on some star formation threshold will re- 
quire fine tuning to keep the threshold at the same location 
with respect to the broad stellar distribution. Models that rely 
on the properties of the gas alone (e.g., Schaye|2004 for gas 
density; Ke nnicutt|1989| for gas stability) will need to explain 
why any accreted gas redistributes itself in such a way as to 
keep the threshold in place for 10 Gyr It is certainly possible 
that the gas responds to the stellar mass already present, as gas 
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flaring and warping seems to indicate ( van der Kruit||2007] l. 
In NGC 4244, the Hi rotation curve drops as the gas density 



begins to decHne, flare, and warp all at about 550" ( |OUing 
[1996 ), somewhat beyond the break radius. A model that de- 
pends on both gas and stellar mass distribution (e.g., Dopita' 
|& Ryder 1994^) might provide a solution, but still would have 
a hard time explaining the sharpness of the break in the very 
old RGB population. Furthermore, all these models may have 
problems explaining why not all galaxies show breaks (e.g., 
NGC 300 has an exponential profile over 10 scale lengths; 
IBland-Hawfliom et al..200 5 ). 

Secular evolution driven by bars and spiral arms can also 
lead to breaks by redistributing angular momentum and hence 
stellar mass (e.g., [DebattistTe t al. 2006). This dynamical re- 
distribution of mass will be the same for all age populations. 
However, if we want to create a break by secular evolution 
of stars alone, the break must be created on a time scale of 
less than 100 Myr. When we account for the fact that the 
gas will also be involved in the secular evolution, we can al- 
low longer time scales, but we then have a similar problem 
as with the threshold models in that star formation must be 
proportional to stellar mass density. Furthermore, the secu- 
lar evolution model requires a significant bar, yet we find no 
evidence for a strong bar in NGC 4244. Our Spitzer 3.6 and 
4.5 micron images (Holwerda et al., in prep.) show an elon- 
gated shape with some twist in the isophotes out to 50" (1/3 
of a disk scale length), but we see no vertical thickening. The 
Hi vel ocity field sho ws some irregularities between 200" and 
300" ( |01Ung||1996| fliat could point to a larger bar, but it is 
not obviously associated with any features in the optical and 
IR images. Finally, the factor 600 decrease in surface density 
at the break (or ^100 compared to the extended inner expo- 
nential disk) makes it unlikely that the break is just a spiral 
arm seen edge-on. This amplitude of arm-interarm contrasts 
has never been observed in face-on galaxies. 

Tidal stripping by galaxy interactions might be another 
route for creating galaxy breaks (and anti-breaks). This mech- 
anism is only effective after multiple interactions, as in a 
galaxy cluster environment ( Gnedin|[2003 ). This is unlikely 



to be the case here as NGC 4244 is the second brightest mem- 
ber of a very loose group with mainly very late-type irregu- 
lar members and NGC 4244 has no detected Hi companions 
within 40 kpc pahlem et al.|2005| ). 

Disk heating and star stripping by pure dark matter sub- 
halos with a ACDM power spectrum is the final model for 
disk breaks we consider here (Kazantzidis et al., in prep.). 
This model has theoretical similarities to the galaxy /cluster in- 
teraction model described above, except that the mass power 



spectrum of interactions is very different. In this model disk 
stars are heated by the constant bombardment of dark matter 
sub-halos whizzing through the disk. The original stars in a 
thin disk will spread into a thicker disk and a diffuse halo de- 
velops. At the disk outskirts, where disk self-gravity is low, 
stars are stripped and a truncation develops. If the stripping 
occurs quickly, independent of stellar age and scale height, 
this model might well explain our observations. 

The stripped stars would form a diffuse, very flattened halo 
around the galaxy with stream-like substructure due to the in- 
teractions with the dark matter sub-halos. This is not unlike 
the structure we see around NGC 4244, with extended compo- 
nents along the major and minor axis of similar surface den- 
sities, which, if present around the whole galaxy, would form 
a very flattened extension around the disk. Unfortunately, the 
number of stars detected in the halo of NGC 4244 is too low to 
measure substructure. Any stellar halo formed through sub- 
halo stripping should not be confused with the stellar halos 
that form from disrupted satellites in hierarchical galaxy for- 
mation. The diffuse halo we detect here is more massive than 
predicted by these kinds of models for galaxies with masses 
similar to that of NGC 4244 ( [Purcell et al.|2007| ). 

It is conceivable that several of the above models are at play 
simultaneously. To have a sharp, deep factor of 600 break in 
the old RGB population one likely needs one of the dynamical 
methods, in which the truncation radius is probably set where 
the stellar disk self-gravity no longer dominates the potential 
and where warps are predicted to start (Saha & Jog 2006). 



Once one of the dynamical methods has created a truncation, 
later infalling gas will respond to the overall potential in such 
a way as to set up a star formation threshold at that location. 
Alternatively, the dynamics of a spiral density wave may in- 
teract with a star formation threshold to set a similar feature 
in the radial profile for young and old stars. Notwithstanding 
this caveat, dynamical break explanations are favored as they 
more naturally explain why profile shapes should be similar 
for both young and old stars. Analysis of similar GHOSTS 
data on a large number of galaxies promises to strengthen the 
constraints on these models. 
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